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ABSTRACT 



Context. Determination of the mass functions of open clusters of different ages allows us to infer the efficiency with which brown 
dwarfs are evaporated from clusters to populate the field. 

Aims. In this paper we present the results of a photometric survey to identify low mass and brown dwarf members of the old open 
cluster Praesepe (age 590*};!! Myr, distance 190+* " pc) from which we estimate its mass function and compare this with that of other 
clusters. 

Methods. We performed an optical (/ c -band) and near-infrared (J and /f s -band) photometric survey of Praesepe covering 3.1 deg 2 . 
With 5cr detection limits of I c = 23.4 and J = 20.0, our survey is predicted to be sensitive to objects with masses from 0.6 to 0.05 M G . 
Results. We photometrically identify 123 cluster member candidates based on dust-free atmospheric models and 27 candidates based 
on dusty atmospheric models. The mass function rises from 0.6 M Q down to 0.1 Mq (a power law fit of the mass function gives 
a=1.8±0.1; £(M)ocM~"), and then turns over at ~0. 1 M Q . This rise agrees with the mass function inferred by previous studies, 
including a survey based on proper motion and photometry. In contrast, the mass function differs significantly from that measured for 
the Hyades, an open cluster with a similar age (r ~ 600 Myr). Possible reasons are that the clusters did not have the same initial mass 
function, or that dynamical evolution (e.g. evaporation of low mass members) has proceeded differently in the two clusters. Although 
different binary fractions could cause the observed (i.e. system) mass functions to differ, there is no evidence for differing binary 
fractions from measurements published in the literature. Of our cluster candidates, six have masses predicted to be equal to or below 
the stellar/substellar boundary at 0.072 M Q . 
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1. Introduction 

Several publications in the past decade have been concerned 
with the mass function (MF) of low mass stellar an d substellar 
popul a tions in open clusters, including <x Orionis (Beja r et al.l 
I2002L ICaballero et all l2007l) the Orion Nebul a Cluster 
([Hillenbrand & Carpenter! 120001 ISlesnicket all \2004, IC 2391 
dBarrado v Navascues et al.l 12004 iBoudreault & Bailer- Jones I 
120091) . the Pleiad es dMoraux et all l2003l lLodieu et al.l | 2007j) . 
and the Hyades dReid & Hawlevl Il999l IBouvier et al.l 12008), 
to name just a few. Studies of relatively old open clusters 
(age > 100 Myr) are important for the following two reasons in 
particular. First, they permit a study of the intrinsic evolution 
of brown dwarfs (BDs), e.g. their luminosity and effective tem- 
perature, which constrains structural and atmospheric models. 
Second, together with younger clusters we can investigate how 
BD populations as a whole evolve and thus probe the efficiency 
with which BDs evaporate from clusters to populate the Galactic 
field. Numerical simulations of cluster evolution have demon- 
strated that the MFs can evolve through dynamical interaction 
dde la Fuente Marcos & de la Fuente Marcosl2000l:lAdams et all 
2002b]). These interactions result in a decrease of the open clus- 
ter BD (and low-mass star) population. This has been observed 



* Based in part on observations carried out at ESO/La Silla, Chile 
under proposal number 078.A-9055(A). 



by IBouvier et al] d2008) from a comparison of the Pleiades 
(120 Myr) and Hyades (625 Myr) mass functions. 

Many earlier studies of the substellar MF have focused on 
young open clusters with ages less than ~ 100 Myr, and in many 
cases much younger (< 10 Myr). This is partly because BDs 
are bright when they are young (lacking a significant nuclear 
energy source, they cool as they age), thus easing detection 
of the least massive objects. However, youth presents difficul- 
ties. First, intra-cluster extinction plagues the determination of 
the intrinsic luminosity function from the measured photom- 
etry. Secon d, at these ages t he BD models have large(r) un- 
certainties dBaraffe et alj|2002l) . Estimates of the substellar MF 
in very young clusters (age <l Myr) might be unre liable due 
to these modelling uncertainties (Cha brier et al . 1 120051) . BDs in 
older clusters suffer less from these problems, but have the dis- 
advantage that much deeper surveys are required to detect them. 

The old open cluster Praesepe is an interesting target con- 
sidering its age and distance. It is located at a distance of 
190^ °pc (based on par allax measurements from the new 
Hipparcos data reduction. Ivan Leeuwenll2009h and has an age 
of 590 + p °Myr (by isochron e fitting in the Hertzsprung-Russell 
diagram; Foss ati et al. I I2008I) . The extinction towards this clus- 
ter is low, E(B - V) = 0.027+0.004 mag dTavlorl 120061) . while 
determinations of the metallicity o f Praesepe yield some dis - 
crepancies: [Fe/Hl = 0.038+0.039 iFriel & Boesgaardl dl992l) : 
+0.13±0.10. lBoesgaard & Budge Q988J); 11+0.03 from spec- 
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troscopy an d 0.20+0.04 from photometry lAn et alj d2007h ; 
+0.27+0. lO. IPace et all d2008l) . lHamblv et all d 1995b presented 
a ~ 19 deg 2 survey of the Praesepe cluster down to masses of 
~O.1M and observed a rise of the MF at the lowest masses. 
They concluded that this implied a large population of BDs. 
A shallow survey complete to 7 = 21.2 mag, R = 22.2 mag over 
800 arcmin 2 uncovered one spectrally confirmed very low-mass 
star or BD (spectral type of M8.5V) with a m odel-dependent 
mass of 0.063-0.084M o dMagazzu et al.lll998l) . A survey over 
the central 1 deg 2 with 10cr limits of # = 21.5, 7 = 20.0 and 
Z = 21.5 mag revealed 19 BD candidates and the first MF de- 
termination of Praese pe down to the subste llar limit, but without 
spectral confirmation dPinfield et al.ll 1997b . Subsequent infrared 
photomet ry of the sample reduced this number to nine can- 
didates (Hodg kin etall Il999l) . lAdams et all (l2002al) presented 
a 100 deg 2 study of Praesepe using 2MASS (Two-Micron All 
Sky Survey) data and Palomar survey photographic plates, from 
which they derived proper motions. They determined the radial 
profile of this cluster but their MF does not reach the substellar 
regime. A more recent proper mo tion survey of Praesepe cov - 
ers a much larger area (300 deg 2 ; |Kr aus & Hillenbrand 2007), 
but does not reach the BD regime either (the limit is ~O.12M ). 
Finally, the most rec ent substellar MF det e rmina tion of Praesepe 
was published by Gonzalez-Garc ia et al. I d2006h and extends to 
a 5<t detection limit of i = 24.5 mag corresponding to 0.050- 
0.055 M . They identified one new substellar candidate, but their 
survey covers only 1 177 arcmin 2 . 

In this paper, we present the results of a program to study, 
in detail, the MF of Praesepe down to the substellar regime . Our 
photometric survey is, as with lGonzalez-Garcfa et al.l d2006l) . the 
deepest so far in optical and near-infrared (NIR) bands, with 5<x 
detection limits of 7 C = 23.4 and J = 20.0 (corresponding to a 
mass limit of about 0.05 M Q ), but covers more than nine times 
the area. Our paper is structured as follows. We first present the 
data set, reduction procedure and calibration in section [2] We 
then discuss our candidate selection procedure in section [3] and 
the survey results in section |4]before discussing the derived MF 
in section [5] We conclude in section [6] 

2. Observations, data reduction, calibration, and 
estimation of masses and effective temperatures 

2.1. Observations 

Our survey consists of 47 Omega 2000 (02k) fields each of 
size 15.4x15.4 arcmin 2 observed in / and K s , plus the same re- 
gion observed in nine 7 C Wide Field Imager (WFI) fields each 
of size 34x33 arcmin 2 . This gives a total coverage of 3.1 deg 2 
observed in all three bands, centred on RA(J2000)=08 /, 40'"04 5 
and DEC(J2000)=+19°40'00". 

The near-infrared (NIR) observations were made on the 
3.5m telescope at Calar Alto, Spain (with observation runs 
of several nights from February 2005 to Januar y 2007). 02k 
(Bail er- Jones et al.l 120001; iBaumeister et al.l 120031) comprises a 
HAWAII-2 detector with 2kx2k pixels over a field of view of 
15.4x15.4 arcmin delivering a pixel scale of 0.45arcsec per 
pixel. The optical observations were carried out with the Wide 
Field I mager (WFI) on th e MPG/ESO 2.2m telescope at La Silla, 
Chile dBaade etalj 19991) during 17-22 March 2007. The WFI is 
a mosaic camera of 4x2 CCDs, each with 2kx4k pixels, covering 
a total field of view of 34x33 arcmin 2 at 0.238 arcsec per pixel. 
All fields were observed in the broad band filter 7 C . A detailed 
list of the fields observed with pointing, filter, exposure time and 
5cr detection limits is given in Table Q] for the NIR data and in 



Table 1. Description of observations with the 02k infrared cam- 
era. 
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Table [2] for the optical data. The passband functions for the fil- 
ters, multiplied with the quantum efficiency of the detectors, are 
shown in Figure Q] 



2.2. Data reduction and astrometry 

The standard data reduction steps (overscan subtraction, trim- 
ming and flat-fielding for the WFI data; dark subtraction and 
flat-fielding for 02k data) were performed on a nightly basis, 
using the ccdred package under IRAH3- For both WFI and 02k 
data we used superflats (obtained by combining science image 



1 The 02k camera suffers from a stray light problem. It appears on 
every image taken with th e camera, fo rming a ring pattern centred in 
the middle of the detector (Nicol 2009). The stray light is removed via 
our global illumination correction and sky subtraction. 
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Table 2. Description of observations with WFI optical camera. 
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Fig. 1. Transmission curve of the filters used in our survey com- 
pared to the synthetic spectrum of a BD with T e ff = 2300 K, log 
g = 4.5 and solar metallicity (NextGen model). The transmission 
curves include the quantum efficiency of the detectors. 



plied an apertu re correction following the technique described in 
iHowelll {1989). An astrometric solution was obtained using the 
IRAF package imcoords and the tasks ccxymatch, ccmap and 
cctran. For each WFI field, this solution was computed for the 
7 c -band image (and for each 02k field using the /-band image) 
using the 2MASS catalogue as a reference. The root mean square 
accuracy of our astrometric solution is 0. 15-0.20 arcsec for both 
WFI and 02k data. For WFI data, the astrometry was performed 
on a CCD-by-CCD basis. 

2.3. Photometric calibration 

To correct for Earth-atmospheric absorption on the photometry, 
we calibrated the infrared data using the J and 7T s -band mag- 
nitudes of 2MASS objects which were observed in our science 
fields. By determining a constant offset between the magnitude 
of 2MASS and our instrumental magnitude, we obtained the zero 
point offset. Since this zero point offset was obtained with ob- 
jects in the same field of view in each science frame, and since 
we found the difference between the 2MASS and 02k pass- 
bands to be insignificant, we did not need to perform an airmass 
or colour correction when reducing our NIR photometry. (That 
is, the determined coefficients were statistically consistent with 
zero.) 

We followed a similar approach for our 7 c -band photom- 
etry, but using observations in our fields for which r and i- 
band magnitudes are available in the Sloan Digital Sky Survey 
(SDSS) catalogue. We first transformed the /-band magnitudes 
of SD SS to 7 c -band mag nitudes using the transformation equa- 
tion of iJordi et alj d2006l) 

I c = i- 0.381 -0.254 X(r-z) (1) 

We then determined the zero point offset between this 7 C mag- 
nitude and our instrumental 7 C magnitude, again using the SDSS 
stars. A further colour correction was not necessary, and as this 
calibration is applied on a field-by-field basis (as with the NIR 
data), an airmass correction was likewise unnecessary. 



frames for each night) for pixel-to-pixel variation correction and 
for correcting the global illumination. For our NIR data, the sky 
background was subtracted using the median-combined images 
for each filter and each field (on a nightly basis). For WFI data, 
we reduced each of the eight CCDs in the mosaic independently 
and in the final step scaled them to a common flux response 
level. We made an initial sky subtraction via a low-order fit to 
the optical data, and for the infrared data by subtracting a median 
combination of all (unregistered) images of the science frames. 
Fringes were visible for the 7 c - band photometry. They were r e- 
moved in the way described by iBailer- Jones & Mundll {200 10. 
Finally, the individual images of a given field were registered and 
median combined. We used the IRAF task daofind to automati- 
cally detect stellar objects in an image by approximating the stel- 
lar point spread function with a Gaussian. We visually inspected 
the images in order to remove from our cluster candidate list any 
extended sources (i.e. galaxies) that were mistakenly identified 
as stars by daofind (see section [33T l. Sources were extracted and 
instrumental magnitudes assigned via aperture photometry with 
the IRAF task wphot. To this aperture photometry we have ap- 



2 A fringe correction frame was created, which is a median combina- 
tion of all science frames in a same filter with the same exposure time. 
This fringe correction frame was scaled by a factor, determined manu- 
ally for each science frame, and subtracted from the science image. 



2.4. Mass and effective temperature estimates based on 
photometry 

After we identify candidates (section 3) we will use the multi- 
band photometry to derive their masses and effecti ve temper- 
atures , T e g. We use the evolutionary tracks from Baraffe et alj 
(1998) and atmosphere models from Ha uschildt et al.l {1999a) 
(assuming a dust-free atmosphere; the NextGen model) to com- 
pute an isochrone for Praesepe for an age of 590 Myr, a dis- 
tance of 190pc, a solar metallicity and assuming zero extinc- 
tion. These models and assumptions provide us with a predic- 
tion of fx, the spectral energy distribution received at the Earth 
(in erg cirT 2 s _1 A -1 ) from the source. We need to convert these 
spectral energy distributions into magnitudes in the filters we 
used. Denoting as Sa(A) the (known) total transmission function 
of filter A (including the CCD quantum efficiency and assum- 
ing telescope and instrumental throughput are flat), then the flux 
measured in the filter is 



£ f A S A (A)dA 
f °°S A (A)dA ' 



(2) 



The corresponding magnitude ma in the Johnson photometric 
system is given by 



niA = -2.5 log f A + ca, 



(3) 
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where ca is a constant (zero point) that remains to be determined 
in order to put the model-predicted magnitude onto the Johnson 
system. We derived this constant for each of the bands 7 C , J and 
K s in the standard way, namely by requiring that the spectrum 
of Vega produce a mag nitude ota of 0.03 in all bands. Using 
the Vega spectrum from IColina et al.l (1 19921) we derive values 
ofc /c = -22.6011, cj = -23.6865 and c Ks = -25.9076mag. 
Applying the two equations above to a whole set of model spec- 
tra produces a theoretical isochrone in colour-magnitude space. 
Note that this procedure only provides us with the "true" mag- 
nitudes of the model spectra, not their instrumental ones. The 
photometric calibration procedure applied to the data converts 
the measured, instrumental magnitudes to the "true" magnitude 
plane where we then compare them with the isochrone. 

Assuming that all our photometric candidates belong to 
Praesepe, we derive masses and effective temperatures from 
these isochrones using our three filter measurements in the fol- 
lowing way. We first normalize the energy distribution of each 
object to the energy distribution of the model using the J filter. 
We then estimate the mass and effective temperature via a least 
squares fit of the measured spectral energy "distribution" (ac- 
tually just two points) to the model spectral energy distribution 
from the isochrone. This involves estimating one parameter from 
two measurements, because mass and T e ff are not independent. 

The above assumption of a dust-free atmosphere is valid 
for T eff > 3000 K, but objects with 3000 K > T eff > 1800K 
are expected to ha ve dust in equilibrium with the gas phase 
dAllard et al.l 1200 1). We therefore perform a second selection 
of candidates (and determination of mass and r e fj) based on 
isochrones p redicted in the same way, but based on evolution- 
ar y tracks of Chabrieret al. (2000) and the AMES-dusty model 
of lAllardet all (120011) . This give us a second dusty model list of 
candidates. A priori some observed stars could appear in both 
lists (and in fact two do), but in our later discussions of the mass 
function we do not mix stars from the two lists but rather make 
separate determinations of the mass function. 

There are various sources of error in the estimation of mass 
and T e ff. These are the photon noise, the photometric calibration, 
the least squares fitting (imperfect model) and the uncertainties 
in the age of and distance to Praesepe. The uncertainties in the 
age and distance are the most significant errors and given rise to 
uncertainties of 0.060+0.010M Q and 1 990+260 K for a substel- 
lar object, 0.072+0.008 M Q and 2 293+201 K for an object at the 
hydrogen burning limit and 1 .000+0.01 7 M Q and 5 300+50 K for 
a solar-type star. 



3. Candidate selection procedure 

The candidate selection procedure for BDs and very low-mass 
stars is as follows (and explained in more detail in the remain- 
der of this section). Candidates were first selected based on 
colour-magnitude diagrams (CMDs) and this further refined us- 
ing colour-colour diagrams. In the third and final selection, we 
used the known distance to Praesepe to reject objects based on a 
discrepancy between the observed magnitude in / and the mag- 
nitude in this band computed with the isochrones and our estima- 
tion of T e ff, To be considered as a cluster member, an object has 
to satisfy all three of these criteria. We make two independent 
selections: one using dust-free and one using dusty atmospheric 
models. 




Fig. 2. Colour-magnitude diagram showing an example of 
the first selection step using the 7 C and J bands. The solid 
lines show the isochrone computed from an evolutionary model 
with a dust-free atmosphere (NextGen model) and the dashed 
lines show our selection band around this. The numbers indi- 
cate the masses (in M ) of objects on the isochrone for var- 
ious 7 C magnitudes. Overplotted are measurem ents from our 
survey of cand idate cluster member s reported in iPinfield et al.l 
(119971) (stars), lAdams et all d2002al) (triangles) [where we in- 
clude objects whic h have a probability of being a real member 
higher than 10%1. iGonzalez-Garcia et al.l d2006l) (squares) and 
iKraus & Hillenbrand! (120071) (circles). 



3. 1 . First candidate selection step: colour-magnitude 
diagrams 

Candidates were first selected from our CMDs by retaining only 
objects which are no more than 0. 14 mag redder or bluer than the 
isochrone in all CMDs. This number accommodates errors in the 
magnitudes, uncertainties in the model isochrones plus uncer- 
tainties in the cluster age and distance estimates. We addition- 
ally include objects brighter than the isochrones by 0.753 mag 
in order to include unresolved binaries. In Figure [2] and [3] we 
show two CMDs where candidates were selected based on 7 C 
vs. I c -J and K s vs. I c -K s . These figures also show low-mass 
cluster member candi dates from previous studies which we de- 
tected in our survey ( Pinfield et al. 1997: lAdams et ai]|2 002a 



iGonzalez-Garcfa et al 



120061: IKraus & Hillenbrand! 120071) . In 
Figure[3] we can observe three structures in this CMD. The two 
structures at 7 C - K s ~ 1 mag and I c -K s ~2 mag are predominantly 
stars (Galactic disk turn-off, and disk late-type and giant stars re- 
spectively) while the structure at 7 C - 7T s ~3mag is mostly com- 
posed of galaxies. From a total of 23 891 objects detected above 
the 5cr detection limit in all filters, 800 are retained as candidate 
cluster members (96.7% are rejected). If we instead use dusty 
model isochrones, then out of the 23 891 objects, 357 are re- 
tained (98.5% are rejected) for our dusty model list. 
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Fig. 3. As Figure[2]but with the I c and K s bands. 



3.2. Second candidate selection step: colour-colour 
diagram 

The second stage of candidate selection involves retaining just 
those objects which lie within 0.24 mag of the isochrone in the 
(single) colour-colour diagram. This value reflects the photomet- 
ric errors plus uncertainty in the age estimation of Praesepe. One 
such colour-colour diagram with the selection limits is shown 
in Figure [4] The two main sources of contamination beside 
field M dwarfs are background red giants and unresolved galax- 
ies (Praesepe is at a Galactic latitude of b-+32.5°). We show 
in Figure [5] the theor etical colours for red gia nts using the at- 
mosphere models of lHauschildt et ail dl999b|) and t heoretical 
colours of six galaxies from Meisenheimer et al. (2009). We see 
that red giants could be a source of contamination in the mass 
range of 0.09-0.2 M and at ~0.7 M , while unresolved galaxies 
should not be a major source of contamination below 0.6 M Q . In 
Figure|5]we see the same three structures as in Figure[3j from top 
to bottom galaxies, disk late-type and giant stars, and Galactic 
disk turn-off stars. Of the 800 objects selected in the first step, 
291 are kept here (63.6% are rejected) assuming a dust-free at- 
mosphere, and 110 out of 357 are kept (69.2% are rejected) when 
using the model for a dusty atmosphere. 



Third candidate selection step: Rejection based on 
observed magnitude vs. predicted magnitude 
discrepancy 



3.3. 



As indicated in section 12.41 our determination of T e ff is based 
on the spectral energy distribution of each object and is inde- 
pendent of the assumed distance. The membership status of an 
object can therefore be assessed by comparing its observed mag- 
nitude in a band with its magnitude predicted from its T e fi and 
Praesepe's isochrone (which assumes a distance). The premise 
is that the predicted magnitude of a background contaminant 
would be lower (brighter) than its observed magnitude and 
higher (fainter) for a foreground contaminant. In order to avoid 
removing unresolved binaries that are real members of the clus- 
ter, we keep all objects with a computed magnitude of up to 
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Fig. 4. Colour-colour diagram used in the second selection step. 
The solid line is the isochrone computed from an evolutionary 
model with a dust-free atmosphere (NextGen model, the masses 
in M for each I c —J colours are shifted up clarity) and the dashed 
lines show our select ion band. Ove r plotte d are th e cluster can- 
didate m embers from Pinfiel det al.1 (1 19971) (stars ), lAdams et al.l 
d2002a|) (triangles), Gonzalez-Gar cfa et al.l d2006l) (squares) and 
from lKraus & Hillenbrand! ^ 2007b (circles), which we detected 
in our survey. 
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Fig. 5. As Figure @] but now showing the theoretical colours 
of six galaxies as thick dotted lines and the theoretical colours 
of red giants as thick solid lines. The six galaxies are two star- 
bursts, one Sab, one Sbc, and two ellipticals of 5.5 and 15 Gyr, 
with redshifts from z=0 to z-2 in steps of 0.25 (evolution not 
considered). We assume that all red giants have a mass of 5 M , 
0.5 < log g < 2.5 and 2000 K < r eff < 6000 K. 
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Fig. 6. Difference between the observed J magnitude and the 
model J magnitude computed from the derived mass and T e g, as 
a function of T e ff. The two vertical lines are at the positions of L0 
and M5 dwarfs (left to right), for orientation purposes. The dot- 
ted line (at -0.753 mag) represents the offset due to the possible 
presence of unresolved binaries, the dot-dash lines represents 
the error on the magnitude determination, and the curved, long- 
dashed lines represent the uncertainties on the age and distance 
of Praesepe. 

0.753 mag brighter than the observed magnitude. We also take 
into account photometric errors and uncertainties in the age and 
distance of Praesepe. This selection procedure is illustrated in 
Figure [6] From 291 objects selected through CMDs and colour- 
colour diagrams in the first two steps, 144 are kept (50.5% are 
rejected) when using the dust-free atmospheres/models, and 35 
out of 110 are kept (68.2% are rejected) when using the dusty 
atmosphere/models. 

After this step, we perform a visual inspection directly on 
the images to reject resolved galaxies and spurious detections. 
This inspection removes 21 and 8 objects from the dust-free and 
dusty selection respectively. 

4. Results of the survey 

We now present the selected candidates, discuss contamination 
by cluster non-members and derive the magnitude and mass 
functions for Praesepe. 



Table 4. Photometric candidates in our survey that are also pho- 
tometric candidates in previous surveys. 



ID 


a 


6 


Alternative name 


Ref/^ 






n 






005 


08 41 08.5 


+ 19 54 02.0 


RIZprl8 


[3] 


010 


08 39 06.9 


+ 19 47 08.0 


J0839069+ 194708 


[1] 








J08390695+ 1947080 


[2] 


011 


08 38 55.46 


+ 19 50 33.3 


J0838554+ 195033 


[1] 








J08385547+ 1950334 


[2] 


012 


08 38 54.19 


+ 19 51 44.6 


J08385420+1951446 


[2] 








J0839030+ 192415 


[1] 








J08390308+ 1924 155 


[2] 


015 


08 39 12.71 


+ 19 30 16.8 


J0839 1272+ 1930 169 


[2] 


016 


08 39 54.39 


+ 19 27 37.1 


J0839544+ 192737 


[1] 








J08395441+1927372 


[2] 


017 


08 39 47.82 


+ 19 28 03.1 


RlZprll 


[3] 


029 


08 42 50.50 


+20 20 03.8 


J0842505+202004 


[1] 








J08425052+2020039 


[2] 


034 


08 42 54.58 


+20 03 36.3 


RIZpr23 


[3] 


035 


08 42 51.96 


+20 05 19.4 


J08425 19+2005 19 


[1] 


038 


08 43 10.76 


+20 01 29.3 


J0843107+200129 


[1] 


040 


08 41 11.05 


+20 22 38.4 


J0841 110+202238 


[1] 


042 


08 40 10.59 


+20 20 50.4 


J0840106+202050 


[1J 








J08401060+2020505 


[2] 


045 


08 39 14.51 


+20 01 19.1 


J0839145+200119 


[1] 


046 


08 39 22.43 


+20 04 54.6 


J0839224+200454 


[1] 








J08392244+2004548 


[2] 


047 


08 38 55.15 


+20 13 08.8 


J0838551+201308 


[1J 








J083855 17+2013089 


[2] 


054 


08 40 53.96 


+20 05 24.3 


J08405397+2005243 


[2] 


062 


08 36 39.46 


+20 22 33.8 


J08363947+2022339 


[2] 


064 


08 36 44.99 


+20 08 45.7 


J08364501+2008459 


[2] 


066 


08 37 11.41 


+20 13 45.8 


J0837 1143+2013459 


[2] 


068 


08 38 08.0 


+20 03 50.1 


J08380800+2003505 


[2] 


070 


08 38 12.44 


+20 08 02.5 


J08381244+2008026 


[2] 


073 


08 38 21.85 


+20 05 35.7 


J08382186+2005356 


[2] 


075 


08 38 39.27 


+19 41 40.1 


J08383929+1941401 


[2] 


081 


08 37 24.48 


+ 1947 11.9 


J08372449+ 1947 120 


[2] 


082 


08 37 02.1 


+ 19 52 07.3 


RIZpr2 


[3] 


101 


08 41 20.32 


+ 18 57 42.9 


J08412034+ 1857430 


[2] 


103 


08 42 19.21 


+19 02 14.8 


J08421923+1902148 


[2] 


108 


08 43 09.0 


+ 1943 11.9 


J08430905+1943119 


[2] 


109 


08 43 01.2 


+ 19 49 59.8 


RIZpr24 


[3] 


117 


08 42 11.47 


+ 19 52 50.2 


RIZpr21 


[3] 


110 


08 43 01.9 


+ 19 54 04.5 


J08430 186+ 1954046 


[2] 


112 


08 42 52.26 


+ 19 51 45.9 


J08425228+ 195 1460 


[2] 


116 


08 42 15.48 


+ 19 48 57.6 


J0842 1550+ 1948576 


[2] 


122 


08 43 08.4 


+ 19 28 06.1 


J08430839+ 1928061 


[2] 


123 


08 43 12.63 


+ 19 34 28.9 


J0843 1265 +1934290 


[2] 



" Objects [1] are from lAdams et ID J2002al). \2\ are from 
iKraus & Hillenbrand! i2007h and [3] are from|pinfield et al] dl997h . 



4. 1 . Selected photometric candidates 

The final selection reveals 123 photometric candidates using an 
isochrone based on dust-free atmospheres, and 27 objects using 
an isochrone assuming dusty atmosphere^. This corresponds to 
~40 and ~9 objects per deg 2 respectively. All our photometric 
candidates are presented in Table [3] Objects are given the no- 
tation PRAESEPE- YYY where YYY is a serial identification 
number (ID). Numbers above 900 indicate candidate members 



3 Two objects in the dust-free atmosphere selection (PRAESEPE- 
089 and -093) were also identified in dusty atmosphere selection 
(PRAESEPE-915 and -917). 



assuming a dusty atmosphere. Only the first 10 rows of the ta- 
bles are shown, all other data are available online. We also note 
in Table [4] which o bjects are candidate c l uster members also 
identifi ed as such by Kraus & Hillenbrand! d2007l) . lAdams et alj 
(I2002al) or lPinfield et al.1 (Il997h . 

Some Praesepe members from previous studies are not de- 
tected in our w ork. This is the case f or the objects from 
iPace et all d2008l) and iFossati et all d2008l) . for example. Since 
those studies foc used on bright obj ects, these stars saturate in our 
science images. dPace et alj r2008 and [Fossati et al J 12008 were 
concerned with chemical abundances of A-type and solar-type 
stars, respectively, while our saturation occurs at ~O.7M .) 
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Table 3. All photometric cluster member candidates of our survey. Table [3] is published in its entirety in the electronic edition of 
Astronomy & Astrophysics. A fraction is shown here for guidance regarding its form and content. 



ID 


a 


6 




J" 


K s a 


M" 


1 eff 


J model 




m 


r") 


[mag] 


[mag] 


[mag] 


[M Q ] 


[K] 


[mag] 


001 


08 40 53.61 


+ 19 40 58.6 


19.19 


16.81 


15.61 


0.089 


2665 


17.00 


002 


08 41 08.8 


+19 43 27.5 


16.14 


14.95 


13.97 


0.219 


3321 


14.86 


003 


08 41 01.6 


+ 19 52 02.5 


16.67 


15.35 


14.38 


0.162 


3189 


15.48 


004 


08 41 12.17 


+ 19 52 48.6 


18.43 


16.39 


15.46 


0.099 


2805 


16.72 


005 


08 41 08.5 


+ 19 54 02.0 


19.02 


16.58 


15.39 


0.088 


2636 


17.06 


006 


08 40 10.74 


+19 40 49.8 


16.97 


15.47 


14.36 


0.132 


3061 


15.95 


007 


08 39 39.56 


+ 19 47 54.3 


17.95 


16.10 


15.07 


0.104 


2860 


16.58 


008 


08 39 43.38 


+19 48 45.7 


16.89 


15.56 


14.68 


0.161 


3186 


15.50 


009 


08 39 55.84 


+ 19 53 14.3 


20.29 


17.50 


16.54 


0.081 


2520 


17.32 


010 


08 39 06.9 


+ 19 47 08.0 


16.51 


15.14 


14.20 


0.155 


3166 


15.57 



" The 1<t uncertainty in the determination of magnitude, effective temperature and mass are the following : Amag = 0.002 mag, AT eff = 140 K 
and AM = 0.1M o for stars (M>0.2M Q ), Amag = 0.01 mag, Ar cff = 230K and AM = 0.05M G for very low-mass stars (0.072 < M < 0.2 M s ), 
Amag = 0.04 mag, AT^g = 420 K and AM = 0.02 M Q for BDs (M < 0.072 M G ). The magnitude 7 mo dei is the predicted magnitude based on photomet- 
ric determination of T eff and mass. 



Not all objects identified by other surveys as brown dwarfs 
or very low mass stellar member candidates - and detected 
in our survey - are membe rs based on our criteria. Th e two 
objects from the work of IGonzalez-Garcfa et al.l (l2006i) . who 
also used photometry in order to select candidate members, 
we detect above our 5<x limit (Prae J084039. 3+192840 and 
Prae J084130.4+190449). Yet both objects are non-members 
based on our selection criteria, because they have I c — J colours 
bluer than our selection band. (Prae J084 130.4+ 190449 is also 
too blue in 7 C - K s for our selection band at I c - K s = 3.0 mag, 
whereas P rae J084039.3+192840 at 7 C - # s =4.0mag lies 
within it.) Gonzal ez-Garcfa et al.l (120061) did not report any 
NIR photometry for these two objects. Although the non- 
membership of Prae J084039.3+ 192840 can be debated (hig h 
membership probability based on IGonzalez-Garcfa et al.l l2006). 
Prae J084130. 4+190449 is most likely an unresolved g alaxy 
(low membership probability; IGonzalez-Garcfa et al.ll2006l) . 

Of the candid ates from the photometric survey of 
iPinfield et al.l (fl997h . seven fall within our survey and are de- 
tected, of which six are identified as candidates by our selec- 
tion cr iteria. The non-selected object is RIZpr6 in Hodgkin et al. 
( 1999). It is bluer than the isochrones in both CMDs in Figure|2] 
and [3] From its positions in the CMDs and in the colour-colour 
diagram in Figure [4] we suspect that this object is an unresolved 
galaxy. 

11 of the the 14 obje cts from a su r vey b ased on proper mo- 
tion and photometry by lAdams et al.l (l2002al) are identified by 
our selection. The objects not recovered fail the observed mag- 
nitude vs. pre dicted magnitude test. On the other hand, 27 cluster 
candidates of Krau s & Hillenb rand (2007) out of 37 detected in 
our survey are selected. T he 10 non-selected objects h ave mem- 
bership probabilities from lKraus & ffillenbrandl d2007t) based on 
proper motion greater than 95%, and are brighter than the 10cr 
detection limit of the publicly available surveys used in their 
work. However, these objects failed our observed magnitude vs. 
predicted magnitude test and some are bluer than our isochrone 
of Praesepe in J - K s . With 7 C - K s colour of ~ 2 mag, we sug- 
gest that these objects are more likely to be disk late-type stars 
or giant stars. 

The 5cr detection limits of our survey are 7 c =23.4mag, 
.7=20. Omag and 7<f s =18.6mag (which correspond to ~0.05M o 
using the dust-free isochrone). However, we cannot expect to 
detect all objects down to these magnitudes. We estimate the 
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Fig. 7. Estimation of the completeness limit for our survey using 
the J band. The solid line is the best linear fit before the turn off, 
the vertical dashed line is the 5 cr detection limit and the vertical 
dotted line is the magnitude at which detector saturation occurs. 



survey completeness by taking the ratio of the number of ob- 
jects detected to the predicted number of detections, the latter 
made by assuming a uniform distribution of stars along the line 
of sight in our survey fields. (This comparison distribution is 
somewhat crude, but it gives an approximate value without mak- 
ing too many assumptions.) The predicted number of detections 
is obtained from the histogram of the number of detections as 
a function of magnitude (Figure |7]i and by observing where the 
distribution drops off compared to a straight line extrapolation. 
Based on this, the completeness of the survey down to the 5<x de- 
tection limit is 90% in 7 C , 88% in J and 87% in K s . The overall 
detection completeness of our survey, from saturation to 5cr de- 
tection corresponding to 0.05 M , is therefore ~87%. In J band, 
we reach a completeness of 95% at J- 19.7 mag, which corre- 
sponds to -0.055 M . 
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Fig. 8. Finding charts of the six new BD candidates of Praesepe 
(/-band). We observed objects very close to PRAESEPE-099 
and -909, although they do not influence the photometry. The 
panels are 50x50 arcsec with North up and East to the left. 

4.2. Substellar candidates in Praesepe 

Six objects in our survey are cluster candidates with theoreti- 
cal masses equal to or below the stellar/substellar boundary at 
0.072 M . We present the finding charts of the six objects in 
Figure [8] In Table [5] we present their coordinates and physi- 
cal parameters. These BD candidates have predicted masses be- 
tween 0.064 and 0.072 M . A spectroscopic follow up (on a 8 m 
class telescope or larger) will be needed in order to confirm or 
refute their membership and their substellar status. 




Fig. 9. / band luminosity function. The solid line histogram 
represents the luminosity function based on a selection using 
a dust-free atmosphere (NextGen model); the thick dotted his- 
togram uses a dusty atmosphere (AMES-Dusty model). The stel- 
lar/substellar limit is at J~ 17.8 mag for both models (0.072 M Q ). 
For reference, the ordinate value of 1 . 1 1 at the histogram peak 
(magnitude bin J = 15.25 mag) corresponds to 13 objects. 

V34 = 5.8 pc than to the 1 1.8 pc error span of the mean cluster 
distance). Therefore, we do not expect contamination by field M 
dwarfs to play a significant role in the determination of the MF. 



4.3. Contamination by non-members 



As mentioned in section 13.21 the two main sources of con- 
tamination are the background red giants, which are the domi- 
nant source at masses of 0.09-0.2 M and ~O.7M , and unre- 
solved galaxies, mostly affecting masses above 0.6 M . Other 
possible contaminants are field M dwarfs and high redshift 
quasars (for instance at z ~ 6; Cab allero et al.l2008T) . However, as 
such quasars have spectral energy distributions similar to mid- 
T dwarfs whereas our faintest candidates are early L dwarfs, 
and given that they are rare (3 .3 quasars at 5.5<z<6.5 in a 
8 deg 2 survey. [Stern et al.ll2007h . the MF should not be affected 
by quasar contamination. 

Let us estimate the contamination by M dwarfs, First, we 
consider that close to the open cluster Praesepe, the space den- 
sity of M dwarfs is uniform. We assume that their density (p) 
drops exponentially with vertical distance from the galactic disk 
(h) such that 



p(h) = p e 



(4) 



assuming a scale height of ho = 500pc. We use the local space 
density (po) for M dwarfs of 5 7-10~ 3 pc~ 3 (from t he Research 
Consortium on Nearby Stars; iHenry et"a"D l2006h . Given the 
Galactic latitude of Praesepe of b-+32. 5 deg and its distance 
of 190pc, the density of M dwarfs near Praesepe should be 
47-10~ 3 pc- 3 . If we define a volume corresponding to the area 
of our survey (3.1 deg 2 , 34 pc 2 ) and use the distance uncertain- 
ties to the cluster (190^'gpc) as its depth, we estimate that we 
have ~19 M dwarf contaminants near the cluster. From a total 
of 150 photometric candidates, we estimate a contamination of 
~13% (or even less, as the cluster depth is presumably closer to 



4.4. Luminosity function and mass function 

We present in Figure|9]the luminosity function of Praesepe using 
the /-band magnitude of the cluster candidates. No correction is 
made for binaries, so this is the system rather than single-star 
luminosity function. 

The mass function (MF), ^(logioM), is generally defined as 
the number of stars per cubic parsec in the logarithmic mass in- 
terval logioM to logioM + c/logioM. Here, we do not compute 
the volume of Praesepe so instead we define the MF as the to- 
tal number of objects in each 0.1 logioM bin per square degree. 
Since we do not make any corrections for binaries we compute 
here a system MF. Our inferred MF is shown in Fig. [10] The 
log-normal form for a MF is 



£(log I0 M) = k-exp 



(log 10 M-log I() M ) 



2o- 2 



(5) 



where A:=0.086, m t 
field system MF 



i o=0.22 Mp and cr=0.57 for the Galactic 
(lChabrierll2003h . Fitting this function to 
both the dust-free and dusty MF data we obtain £=5.9+3.1, 
m o =O.15±O.O5M and cr=0.51+ 12. Figure [101 shows this re- 
sult. If we instead fit a power law (ISalpeterl (1955)) 



f(log 10 M) = *.M- 



(6) 



from the highest mass bin to the turn over at 0. 1 M , we obtain 
or=1.3±0.2 (corresponding to £(M) oc M~ 2 3 ). If we exclude the 
two bins between 0.1 and O.15M (possible contamination by 
red giants) and the two highest bins (possible incompleteness), 
the fit gives a=0.8±0.1. 
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ID 


a 


8 




J 


K s 


M 


Tefi 


J model 


Spill 






C") 


[mag] 


[mag] 


[mag] 


[Mo] 


[K] 


[mag] 




055 


08 41 04.5 


+20 14 58.0 


21.61 


18.29 


17.12 


0.068 


2250 


17.97 


L0 


096 


08 41 13.48 


+18 59 05.1 


21.06 


17.85 


16.82 


0.072 


2335 


17.75 


M9 


099 


08 41 45.16 


+ 19 18 07.7 


21.30 


17.98 


17.01 


0.068 


2249 


17.98 


L0 


909 


08 39 29.94 


+20 11 40.3 


20.11 


17.63 


16.67 


0.069 


2259 


17.95 


LO 


910 


08 40 34.00 


+20 14 56.2 


20.08 


17.60 


16.65 


0.069 


2261 


17.94 


LO 


915 


08 38 51.77 


+ 19 00 21.6 


20.28 


17.67 


16.68 


0.068 


2238 


18.01 


LO 



" Spectral type expected based on T e g and the temperature scale of iKraus & Hillenbrand i2007h . with LI set to 2100 K. The error on this 
estimation is one subclass. 
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Fig. 10. Mass function based on our survey photometry. Points 
with error bars represent the MF based on a selection and mass 
calibration assuming a dust-free atmosphere, whereas the open 
circles with error bars are the MF based on the dusty atmosphere 
model. We also overplot the log-normal and the power law MF 
fitted to our data (both solid line). Error bars are Poissonian aris- 
ing from the number of objects observed in each bin. The vertical 
thin dotted lines are the mass limits at which detector saturation 
occurs in the I c , J and /T s -bands (from left to right). The vertical 
thick dashed line is the mass at the 5cr detection limit (complete- 
ness of ~87%). For reference, the ordinate value of 0.932 at the 
histogram peak (mass logioM=-0.85 [O.142M ]) corresponds 
to 27 objects. The two dusty data points have been shifted to the 
right by logioM=0.05 for clarity. 



5. Analysis and discussion of the stellar and 
substellar mass function of Praesepe 

Our MF of Praesepe (Figure |T0b shows a rise in the number of 
objects from 0.6 M down to 0.1 M , and then a turn-over at 
~0.1 M . This turn-over is not due to incompleteness since it oc- 
curs well above the 5<x detection limit corresponding to 0.05 M . 
This behaviour is confirmed by the luminosity function in Figure 
[9] which shows a rise from 7=13 to 16 mag (with candidates ob- 
tained using a dust-free atmosphere) and a drop at J- 17 mag 



(seen with both types of candidates). To help the analysis of 
these features in the mass function, we compare in Figure[TT1the 
mass functions of Praesepe obtained from several studies plus 
the MF for the old open cluster Hyades (age of 625 Myr). 

The rise in our MF of Praesepe is also present in the MFs ob - 
tained in the three previous studies of Baker & Jameson (2009), 
IKraus & Hillenbrand! d2007l) and lHamblv et all d 19951). On the 
other h and, we do not see this rise in the MF of lAdams et alj 
(2002a). However, their MF is based on objects with a mem- 
bership probability higher than only 1 % and within a radius of 
3.8 deg. Due to use of such a low probability threshold for se- 
lection, we expect that most of the objects used in the MF de- 
termination are si mply field stars (wh ich is their own conclu- 
sion in section 5.4; Ada ms et al.l2002a|). so further compa r ison is 
not warranted. As for th e MFs of [Gonzalez ^Garcfa et ail |2006) 
and lPinfield etail dl997l) . since the highest mass bins are ~0. 1 1 
and ~O.15M (respectively), the rise observed from 0.6 M to 
0. 1 M cannot be discussed. 

Only fou r MFs, in addition t o our work, reach masses be- 
low .1 M : iBaker & Jameson! (120091). iGonzalez-Garcfa et alj 
d2006b . IPinfield et al l (fl997h a nd lHamblv et all (fl995l) While 
the MFs of IBaker & Jameson] (120091) and lHamblv et al] (1 19951) 
show a turn-over at 0.1 M , the one obtained by Pinfield et al. 
( 1997) does not. On the contrary, it presents a sudden rise at the 
stellar/substellar limit (with a ratio of ~5 in the number of ob- 
jects at the mass bin at 0.07 to the number in the bin at 0. 1 1 M ). 
They used RIZ photometry for their survey, but not all objects 
were observed in all bands, resulting in just one colour avail- 
able f or membership determination in some cases (Pinfield et al. 
1 19971) . From an analysis of MFs of other clusters and using 
a mul ti-band photometric survey, iBoudreault & Bailer- Jones I 
(2009) have shown that use of a narrow spectral coverage with 
few filters can lead to high contamination by field M dwarfs, 
and thus an apparent rise in the MF in the low mass regime. We 
suggest that this is the r eason for the a ppare nt rise at the low- 
mass end of the MF in IPinfield et al.l (1 19971) (who also noted 
that only one colour is availa ble for many objects in their two 
lowest bins). As for the MF of Gonzal ez-Garcfa et al] d2006l) . as 
they only have three points we cannot comment on any possible 
trend. 

Although there are some discrepancies between the dif- 
ferent MFs of Praesepe from previous works and our MF, 
none agre es with the MF of the Hyades (~625 Myr) ob- 
tained by iBouvier et ail (|2008fl in which the MF is ob- 
served to turn-over and decrease already at 0.35 M . This 
is surprising, since Praesepe and the Hyades share a com- 



4 Like the MF of Praesepe we present, the MF of the Hyades pre- 
sented bv lBouvier et alj d2008h is a system MF (no correction for bina- 
ries). 
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which have no data here. 



parable age, size and mass: they have ages of 590|^°My r 
dFossafi et al] 120081) and 625±50Myr dBouvieret all 12 008), 
tidal radii of 1 1.5+0.3 pc ( 3.5+0.1 deg. iRraus & Hillenbrand! 
|2007|) and 10.3 pc (12.5 deg, iBouvier et al.ll2008b . and masses 
of 550+40 Mr ((Kraus & Hillenbra nd 2007) and about 400 M 
dBouvier et all 120081) . respectively. Therefore, we can expect 
that the potential well is the same (at least today). Only the 
metallicity may be slightly different, assuming the most re- 
cent measurement for Praesepe: [Fe/H] = +0.27+0. 10 for th e lat- 
est metallicity measurement of Praesepe (Pace et al. 2008) and 
[Fe/H] = +0.14+0.05 for the Hyades Bouvier et all 120081). al- 
thoug h a metallicity as low as +0.038+0.039 dFriel & B oesgaard 
1992) has been reported for Praesepe. It is unclear how this 
metallicity difference could explain the significantly different 
mass functions. 

It is a priori possible that different binary mass fractions in 
Praesepe and the Hyades could account for the difference in their 
observed (i.e. system, rather than star) mass functions. The bi- 
nary fract ion in Praesepe for different mass intervals was ob- 
tained bv lPinfieldetaiTd2003l) : 17^% for 1. 0-0.6 M G , 31^% 
for 0.6-0.35M o , 44+6% fo r 0.35-0.2 M Q and 47 + }^% for 0.1 1- 
0.09 Mq. As for the Hyades. iGizis & Reidldl995l) observed a bi- 
nary fraction of 27+16% for their sample of stars (<0.4 ), which 



is consistent with anothe r determination of the Hyades binary 
fraction of 30±6% from iPatience et ail (H998) (for a primary 
mass of ~0.6-2.8Mq). From these figures we see no significant 
difference in the binary fractions of the two clusters (even if pri- 
marily because the uncertainties are quite large), so this cannot 
be used to explain the difference between in their mass functions. 
Of course, if the typical mass ratio in a binary system is differ- 
ent in the two clusters then this may be able to account for some 
difference in the mass functions, but their is also no evidence to 
support (or refute) this. 

A distinction between the two cluste rs could be the spa- 
tial distribution of the members. Indeed, Holland et al. (2000) 
observed that the Praesepe cluster might be composed of 
two merged clusters with different ages, one main cluster of 
630 M and a second subcluster of 30 M G . It was even proposed 
that faint low-mass members of the s ubcluster could appea r 
as Praese pe brown dwarf ca ndidates (Chappe lle et al.l 12005). 
However, Ada ms et al.l d2002al) did not find evidence of a sub- 
cluster in Praesepe. Based on the spatial dis tribution of the main 
cluster and subcluster from Holl and et all d2000l) . our survey 
only overlaps the main cluster. In addition, a collision between 
two clusters could not explain alone an increase of the MF down 
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to 0. 1 M Q , as such a collision would rather remove low-mass 
member of the clusters. 

By comparing the MF of the Hyades with the one of the 
Pleiades (~ 1 20 Mvr). iBouvier etaT] (120081) concluded that dy- 
namical evolution was responsible for the deficiency observed in 
the very-low mass star and BD regime in the Hyades. However, 
this deficiency is not seen in Praesepe. One possible implica- 
tion is that Praesepe has been less affected by dynamical evolu- 
tion, i.e. evaporation of low mass members which are expected 
to have higher speeds based on equipartition of energy. On the 
other hand, if dynamical evolution has affected Praesepe in the 
same way, then it cannot have had the same initial mass function 
and/or initial conditions as the Hyades. Dynamical interaction 
between one of these clusters and another object (such as an- 
other open cluster in the past) could explain the discrepancies 
between the two MFs. 

6. Conclusions 

We have presented the results of a survey to study the mass func- 
tion of the old open cluster Praesepe. The survey consisted of 
optical 7 c -band photometry and NIR J and ^ s -band photometry 
with a total coverage of 3.1 deg 2 , down to the substellar regime, 
with a 5cr detection limit corresponding to 0.05 M (the detec- 
tion completeness to this level is ~87%). 

Our final sample yields 123 photometric cluster member can- 
didates based on a selection assuming a dust-free atmosphere 
and 27 photometric cluster candidates based on a selection as- 
suming a dusty atmosphere. We estimate the contamination by 
field M-dwarfs to be 13% or less. Among our cluster candidates, 
six objects have theoretical masses equal to or less than the stel- 
lar/substellar boundary at 0.072 M . 

We observed that the MF of Praesepe is characterized 
by a rise in the number of objects from 0.6 M Q down to 
0.1 M , followed by a turn-over in the MF at -0.1 M Q . The 
rise is in agreement with the Praesepe MFs derived in sev- 
eral previous studies dHamblv et al.lll995HKraus & Hillenbrand 
120071; iBaker & Jameson 2009) but disagrees with lAdams et al. 
d2002al) . 

We have compared the mass function of Praesepe with one 
derived for the Hyades and have observed a significant differ- 
ence: while the Hyades has a maximum at 0.35 M , Praesepe 
has a maximum at a much lower mass, 0.1 M . Assuming that 
they have similar ages (as main sequence fitting suggests), we 
conclude that the clusters either had different initial mass func- 
tions or that dynamical interaction has modified the evolution 
of one or both. More specifically, in the latter case, dynamical 
evaporation does not seem to have influenced the Hyades and 
Praesepe in the same way. A difference in the binary fraction or 
mass ratios could also cause a difference in the mass functions, 
but determinations of these are not yet precise enough to suggest 
any difference. 
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